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Abstract. We consider a nonlocal aggregation equation with non- 
linear diffusion which arises from the study of biological aggrega- 
tion dynamics. As a degenerate parabolic problem, we prove the 
wcll-posedncss, continuation criteria and smoothness of local solu- 
tions. For compactly supported nonnegative smooth initial data we 
prove that the gradient of the solution develops L!^-norm blowup 
in finite time. 



1. Introduction and main results 

In this paper we consider the following evolution equation with non- 
linear diffusion: 

dtu + d^{ud^K *u) = rd^iu^d^u), {t, x) e (0, oo) x R, 
m(0,x) = no(x), 

where K is an even function of x and has a Lipschitz point at the 
origin, e.g. K{x) = e~'^'. Here * denotes the usual spatial convolution 
on M. The parameter r > measures the strength of the nonlinear 
diffusion term. The function u = u{t,x) represents population density 
in biology or particle density in material science, (see, for example, 
[56][26][30][3T][16][57][^[M][2Z][1Z] [SI]). This equation was 
used in the study of biological aggregation such as insect swarms, fish 
schools and bacterial colonies. It is first derived by Bertozzi, Lewis 
and Topaz as a modification (more precisely, the addition of the 
density-dependent term on the RHS of fll.ip ) of an earlier classical 
model of Kawasaki According to [SB], this modification gives 

rise to biologically meaningful clumping solutions (i.e. densities with 
compact support and sharp edges). For other similar one-dimensional 
models and their biological applications we refer the readers to [16] [26] 
[60] [55] [29] [32] [33] [31] [M] [19] [50] and the references therein for 
more extensive background and reviews. 

To understand the biological meaning of each term, one can rewrite 
(II. ip as the classical continuity equation 

dtu + d,{uv)=0, (1.2) 
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where the velocity v is related to the density u by 

V = dx{K * u) — rudxU 

Here Va is called the attractive velocity since as explained in [SB] in- 
dividuals aggregate by climbing gradients of the sensing function s = 
K * u. Due to the spatial convolution Va is a nonlocal transformation 
of the density u. The second term Vd is called the dispersal (repulsive) 
velocity and is a spatially local function of both the population and the 
population gradient. Biologically Vd represents the anti-crowding mech- 
anism which operates in the opposite direction of population gradient 
and decreases as the population density drops. These constitutive rela- 
tions of fa, Vdi u are natural in view of the basic biological assumption 
that aggregation occurs on a longer spatial scale than dispersal. 

In the mathematics literature, the aggregation equations which have 
similar forms to f 1 1.1 1) have been studied extensively [TT] [T2] [13] [E] 
[IS] [3H] [57]. The case of (11.11) with r = and general choices of 
the kernel K was considered by Bodnar and Velazquez [13]. There 
by an ODE argument the authors proved the local well-posedness of 
(II. ip without the density- dependent term for initial data. For a 
generic class of choices of the kernel K and initial data, they proved 
by comparing with a Burgers-like dynamics, the finite time blowup of 
the L!^-norm of the solution. Burger and Di Francesco [H] studied a 
class of one-dimensional aggregation equations of the form 

dtu = dx {udx{a{u) — K *u + V)) , in (0, oo) x R, 

where V : R ^ M is a given external potential and the nonlinear 
diffusion term a{p) is assumed to be either or a strictly increasing 
function of p. In the case of no diffusion (a = 0) they proved the 
existence of stationary solutions and investigated the weak convergence 
of solutions toward the steady state. In the case of sufficiently small 
diffusion (a(p) = ep^) they proved the existence of stationary solutions 
with small support. Burger, Capasso and Morale [I5] studied the well- 
posedness of an equation similar to (II. ip but with a different diffusion 
term: 

dtU + V ■ {uVK *u) = dzv{uVu), in (0, T) x R"^. 

For initial data Uq E L^(R'^) n L^(R'^) with e HKR"^), they proved 
the existence of a weak solution by using the standard Schauder's 
method. Moreover the uniqueness of entropy solutions was also proved 
there. In connection with the problem we study here, Laurent [3H] has 
studied in detail the case of (II. ip without density- dependent diffusion 
(i.e. r = ) and proved local and global existence results for a class 
of kernels K with H^{s > 1) initial data. More recently Bertozzi and 
Laurent [11] have obtained finite-time blowup of solutions for the case 
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of fll.ip without diffusion (i.e. r = 0) in M.'^{d > 2) assuming compactly 
supported radial initial data with highly localized support. Li and Ro- 
drigo [H] [12] [13] studied the case of fll.ll) with fractional dissipation 
and proved finite-time blowup or global wellposedness in various sit- 
uations. We refer to [H] [l5] for the cases with singular kernels and 
further detailed studies concerning sharp asymptotics and regularity of 
solutions. We also mention that Bertozzi and Brandman [7] recently 
constructed fl weak solutions to fll.ll) in [d > 2) and with no 
dissipation (r = 0) by following Yudovich's work on incompressible Eu- 
ler equations [HI]. We refer the interested readers to [55] [22] [32] [SS] 
[SI] [SS] [IH] [SO] and the references therein for some further rigorous 
studies. 

From the analysis point of view, equation fll.ip is also connected 
with a general class of degenerate parabolic equations known as porous 
medium equations, which takes the form 

dtu = {t,x) e (0, oo) X M, (1.3) 

where m is a real number. These equations describe the ideal gas flow 
through a homogeneous porous medium and other physical phenomena 
in gas dynamics and plasma physics [IE] [62] [5]. Oleinik, Kalashnikov 
and Chzhou [52] proved the existence and uniqueness of a global weak 
solution of fll.Sp . Holder continuity of the weak solution is established 
in [1] [37] [2H] [20]- As for higher regularities, the Lipschitz continu- 
ity of the pressure variable v = ^^^^u is shown in [1] [2S] [ID] [2] [21]- 
Aronson and Vazquez [6] showed that v becomes C°° outside of the 
free boundary after some waiting time. For derivative estimates of the 
solution u, a local solution of 01.30 in W^'°°(M.) is constructed by Otani 
and Sugiyama [SS]- In the case of fll.3p with m being an even natu- 
ral number, Otani and Sugiyama [51] proved the existence of smooth 
solutions. We refer the interested readers to P [3] [9] [8] [IE] [IT] [21] 
|23] [21] [IE] [in] [22] [36] [59] and references therein for more detailed 
studies and expositions. 

Our starting point of the analysis of equation fll.ll) is to treat it as a 
degenerate parabolic problem with a nonlocal flux term. We focus on 
constructing and analyzing classical solutions of fll.ll) . The bulk of this 
paper is devoted to proving the existence and uniqueness of classical 
solutions to fll.ip . which is done in section 2 and part of section 3. 
In the final part of section 3, we prove the continuation and blowup 
criteria of solutions. In the last section we prove that any smooth initial 
data with compact support will lead to blowup of the gradient in finite 
time. The analysis developed in this work can be extended to treat the 
d-dimensional {d > 2) case of fll.ll) which we will address in a future 
publication. 
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We now state more precisely our main results. The first theorem 
establishes the existence of smooth local solutions for initial data which 
is not necessarily nonnegative. 

Theorem 1.1 (Existence and uniqueness of smooth local solution). 
Assume Uq G f]^^gif™(M) (uq is not necessarily nonnegative). Then 
there exists a positive Tq = To(||mo||2 + ||<9xiio||oo) such that fll.ll) has a 
unique solution u G C°°{[0,To] x R). In particular u G C([0,To] fl H^) 
for any k >0. If uq > 0, then u{t) > for all < t < Tq. 

Remark 1.2. The assumptions on the initial data Uq in Theorem 11.11 
can be weakened significantly (see for example Theorem 13. 2p . However 
in order to simplify the presentation, we do not state our theorems here 
in its most general form. 

Our second theorem gives the blowup or continuation criteria of so- 
lutions. Roughly speaking, it says that all the L^-norm of the solution 
cannot blow up and we can continue the solution as long as we have a 
control of the gradient of the solution. 

Theorem 1.3 (Blowup alternative and growth estimate of L^-norm). 

Assume u G C°°([0,T) xM) is a maximal- lifespan solution obtained in 
Theorem \l.l[ Then either T = +oo or T < +oo and 

lira \\dxu{t)\\oo = +oo. 

For any 2 < p < oo, there exists a generic constant C such that 

lk(t)||p< ll^ollpe^^*, VtG [0,T). 

// in addition uq > 0, then u(t) > for all t G [0, T), and we also have 
the p-independent estimate for all 2 < p < +oo.- 

\\uit)\\p < llMollpCxp {CtWuohe^') , Vt G [0,T). 

In particular if p = +oo, then 

\Ht)\\oc< ||no||ooexp(Ct||no||2e^*), VtG [0,T). 

The next theorem states that if we assume the initial data has a little 
bit more integrability, then the local solution will inherit this property. 
Note in particular that the L;^-norm of the solution is preserved for all 
time if the initial data is nonnegative and in L^(M). 

Theorem 1.4 (Local solution in Lp{R)). Assume uq G f]'^^QH^(R) 
and Uq G Lix^) /^'^ some 1 < p < 2. Then the local solution obtained 
in Theorem also satisfies u G C([0, T], L^.). // in addition Uq > 0, 
then u{t) > for any t G [0,T]. // also p = 1, then \\u{t)\\i = ||uo||i, 
i.e. L\.-norm of the solution is preserved. 

The last theorem states that any solution with smooth nonnegative 
initial data will blow up in finite time. 
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Theorem 1.5 (Finite time blowup). Let Uq G C^(]R) and Uq > 0. 
Assume Uq is not identically 0. Then there exists time T < oo such 
that the corresponding solution with Uq as initial data blows up at time 
T in the sense that 

lim ||9:i.M(t)||oo = +00. (1.4) 

However all the L^-norm of u remain finite at the time of blowup, i.e.: 
sup ||M(t)||p < C(T) < cx),V 1 < p < oo. (1.5) 

0<t<T 

Remark 1.6. In the case of (11. ip without diffusion (i.e. r = in (11.11) ). 
if we take the kernel K{x) = e~'^', then the result of Bodnar and 
Velazquez [13] says that any solution with initial data Uq satisfying a 
slope condition will blow up in finite time in the sense that ||M(t)||oo 
blows up. This is highly in contrast with our result here when the 
diffusion term does not vanish. In this case the solution will blow up 
at the level of the gradient, i.e. we have ||Vu(t)||oo tends to infinity 
while all the other L^-norm remain finite as t — > T, where T is the 
blowup time. 

Notations. Throughout the paper we denote = I>^.(M) for 1 < 
p < oo as the usual Lebesgue space on M. We also write || ■ ||p = || ■ Hig- 
For s > 0, s being an integer and 1 < < oo, W^'^ = PV^'^(M) denotes 
the usual Sobolev space 

W:'^ = {/ G S'iR) : ll/llvy.P = W^'JWp < 

0<j<s 

When p = 2, we denote = H^{R) = W^'P{M.) and || ■ as its 
norm. Occasionally we shall use the Sobolev space of fractional power 
if|(M) whose norm can be defined via Fourier transform: 

ii/ik» = ii(i + iei)7(0iiL|. 

For any two quantities X and Y, we use X < Y or Y > X when- 
ever X < CY for some constant C > 0. A constant C with sub- 
scripts implies the dependence on these parameters. We write A = 
A{Bi, ■ ■ ■ , Bk) when we want to stress that a quantity A depends on 
the quantities -Bi, ■ ■ ■ , Bk. 

From now on we assume r = 1 in (II. ip without loss of generality. 
Same results hold for any r > 0. 

2. The regularized equation and its wellposedness 

Since (II. ip is a degenerate parabolic equation, in order to construct 
a local solution, we have to regularize the equation. To this end, we 
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consider the following regularized version of fll.ll) 

dtu = (u^ + e)d^^u + 2u{d^uY - d^{uK * d^u), 
u(0, x) = Uq{x). 

Here e > is a parameter. We are going to prove the following 

Proposition 2.1 (Local solution of the regularized equation). Assume 
Uq G f]'^^Q H^(M) (uq is not necessarily nonnegative). Then there 
exists a positive Tq = Todl^^jWoHoo + ll'^olh) (Tq is independent of e) 
such that fl2.ip has a unique solution u'^ G C°°([0,To] x M) for any 
e > 0. In particular u'^ G C([0,To]; i/^) for any k >0. 

Corollary 2.2. Assume Uq G f^^^Q i/^(M) and uq G L^(^) for some 
1 < p < 2. Then the local solution obtained in Proposition \2.1\ also 
satisfiesu'' G C([0, Tq], L^). AndTo can be chosen to beT^ = To(||mo||p+ 

ll^x-Molloo)- 

In the case uq G fXn=o with uq > 0. The regularized equation 

admits a global solution. 

Proposition 2.3. Assume Uq G f]^^Qif^(M) with uq > 0. Then for 
anyT > (12.11) has a unique solution u"^ G C°°([0, T] xM). In particular 
u" G C([0,T]; H^) for any k > 0. Also u'{t) > for any t>0. If m 
addition uq G Ll{R) for some 1 <p <2, then u" G C([0,T]; L^). 

2.1. Proof of Proposition [27TI Our proof of Proposition 12. II is rem- 
iniscent of the energy method used by Otani and Sugiyama [M] 
where they dealt with the porous medium equation (II. 3p . Denote the 
set 

5^ = {t; G C([0,r]; Hf+\m)) : d^^v^dtV G L\[Q,T]- i/f (K))}- 

As a very first step, we shall show the local existence of the solution to 
(12. ip in B^. At this point, we need the following lemma from |54j . 

Lemma 2.4. Consider the initial value problem 

dtU = [u^ + e)d^^u + h{t, x), 

n(0) = Uq, 

where h G L^(0,T; H^''), Uq G H^''~^^{M.), there exists a unique function 

u = u{t, x) G which satisfies, for some constant Ci = Ci(||Mo||j:^2fc+i, 2,2^1^2*!, k, e), 

that 

LOO jj^k+i -\- e^\\u\\^^2k+2 < Cl- 

For any hi, h2 G Kj^ = |i; : ||^'||L2(o,T;_H"2fe) < R} , there exists a 
constant C2 = C2{R, k, e) such that 

\\Ul — M2II roo rr2fe + l + £2 \\Ul — M2II r2rr2fc + 2 < C2c'^^'^||^l ~ h2\\ ffik . 

t X fx t X 

Proof. See [51]. □ 
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Based on this lemma, we establish the local existence of solutions of 
the regularized equation. Note however that the time of existence of 
the solution depends on the regularization parameter e. We have the 
following 

Lemma 2.5 (Local existence of the regularized equation). Let Uq G 
H^^+^(m), k>2. There exists Tq = To(||Mo||//2fe+i, ^, e) > 0, such that 
(12. ip has a unique solution u G B^^. 

Proof. Introduce the auxiliary equation 

dtu = {u^ + e)d^^u + h{t, x) 

U{0) = Uq 

By Lemma [Ml for any Uq G iff+^M), h G L'^{<d,T; Hf), Q has a 
unique solution G B^. Next define the map 

It suffices to show that for some suitable R and Tq, is a contraction 
from the set = {v e L'^{0,To; Hf); ||i;||L2(o,To;//|fc) < R] into 

itself. First we show that </> maps into . Take any h G -fC^, by 
the fact that H^{^) is an algebra when s > |, we have 

\\(j){h)\\H2k < \\uh\\H^k\\uh\\]j2k+i + \\d^K * dxUhUhWnik + \\K * d^Uhd^UhWH^k 
<C{\\ 

where C is a generic constant. Then Lemma [2.41 implies that 

l|0(^)l|L2(O,T;H2fc) < C'^^(ll^/i|lLJ«(0,T;H2fc+i) + II II Lf> (0,T;/f2'=+i) ) 

< T^C3(||Mo||^2fc+i,/?, /c,e). 

This shows that if we take T < C|/i?^, then maps K'^ into K^. 
Similarly we can show that is Lipschitz, i.e. for any hi, /i2 G K^, we 
have 

||0(/ll) - 0(/l2)||L2(O,T;H2fe) < C4r5e^^^||/li - /l2 1| L2(o,T; //2.) , 

where C4 = C4(||no||j^2fc+i, i?, /c, e). Therefore if T is sufficiently small, 
then is a contraction on Kj^. The lemma is proved. □ 

Lemma [2.51 is not satisfactory since the time of existence of the solu- 
tion depends on the regularization parameter e. The following lemma 
removes this dependence, and at the same time, weakens the depen- 
dence on the initial norm. 

Lemma 2.6. Let uq G iJ^fe+i^]^)^ k > 2. Then there exists Tq = 
To{\\dxUo\\oo+\\uQ\\2) such that for any e > 0, (12. ip has a unique solution 
G //no > 0, then u'{t) >0 for anyO<t< Tq. 
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Proof. By Lemma 12.51 we can continue the solution as long as we can 
control the H^^'^^-norm of u'^. In the following we give the a priori 
control of u'^. As we will see, the H^^~^^-noTm of will stay bounded 
on a time interval [0,Tq] for a certain small Tq depending only on the 
norms ||9a;Uo||oo, ll^olh- We are then left with the task of estimating the 
various Sobolev norms which we will do in several steps. For simplicity 
of notations we shall write u"^ as u throughout this proof. 

Step 1: L^-norm estimate for 2 < p < +oo. Multiply both sides of 
(12. ip by \u\P~'^u and integrate over x, we have 



p dt 



= —e{p — 1) J \u\^ '^uldx — {p — 1) J ul\ufdx 

+ {p-l) j {dxK * u)u\u\P-^dxUdx (2.4) 

<-(P-i)/i"r(-.-^) 



2 

dx 



+ J \u\P-\d,K*uydx 

<^\\dxK*u\\i\\ur-' 

< ^\\dxK\\l\\u\\i. 



Gronwall's inequality implies that 

\Ht)\\p < hoWp ■ exp (^^^\\dxK\\l?j , V2 < p < +CX). (2.5) 
In particular we have 



\Ht)\\2 < llnolhexp \^^\\dxK\\ltj . (2.6) 

To control L!^-norm, we are going to choose Tq < -^^qj^^- This implies 
that ||tt(t)||2 has a uniform bound on any such [0,To]. Now use again 
the RHS of (EJ) to get 

^ ^ "u{t)\\l <ip-l) J {dj< * u)u\uf-''dxudx 



p dt 

p — 1 

p 

p — 1 
p 



idj{ * u)dx{\uf)dx 
{dxxK * u)\u\^dx. 
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Now we use the fact that K{x) = e^'^'' and therefore dxxK{x) = 
—25{x) + e"!"^!, and this gives 



< ^- / u\u\Pdx + / (e^'^l * u)\u\Pdx 



pdt ^ p J p 

< 1)|| IIP ^~^\\ IIP 

p ^ p ^ 



This imphes 

d 2ip — 1) P ~ 1 II / \ II 

37 F(^) P ^ F P F oo H \\u[t)\\p. 

dt p P 

Integrating over [0,t], we obtain 

\Ht)\\p < \\uo\\p + C [ {l + \\u{s)\\^)\\u{s)\\pds. 
Jo 

Letting p ^ oo, we obtain 

||^i(i)||oo < ||%||oo + C (1 + \\u{s)\\oo)\\u{s)\\oods. 



Now it's easy to see that if we choose Tq sufficiently small depending 
on ||mo||2, ll^iolloo, then we have 

sup \\u{t)\\r < C\\uo\\r, V2<r<oo. (2.7) 

0<t<To 

Step 2: Control of II I loo- At this point, one can appeal directly 

to the maximum principle (see for example Theorem 11.16 of [M])- We 
give a rather direct estimate here without using the maximum principle. 
Differentiating both sides of fl2.ll) w.r.t x and denoting v = u^, we have 

dtv = evxx + u^Vxx + Quvvx + + G, (2.8) 

where 

G = dxx{uK * v). 

Multiplying both sides of (12.81) by |t;|^~^w and integrating over x, we 
obtain 

- — ||f(t)||p = j u'^Vxx\v\^~'^vdx + e j Vxx\v\^~'^vdx + Q j uvx\v\^dx 
+ 2 [ \v\P-^^dx+ [ G{x)\v\f-\dx. (2.9) 



10 DONG LI AND XIAOYI ZHANG 

Integrate by parts and we have 





'"^vdx = 






UVx\ 


v\^dx = 


1 


-f 


1 + 


P J 




~'^vdx = 


-e{p 


-1) 



Plugging the above estimates into (12.91) . we get 

^ d ,, , , , 

Now we compute 



.u,„s<2/i„r.../GWK,rw 



f G\v\''~\dx= I d^^{K * vu)\v\''~\dx 

< j \u\\v\^dx + j {e ^''^*v)u\v\^ "^vdx + \\v\\l\\v\\oo 

^ll^'llp(hl|oo + IK'lloo). 



Finally we obtain 

^|lK'WII^<IK'lC2+IK'll^(IIHIoo+||^^||oo) 
^ ll^llpdl^IlL + ll^lloo + Iklloo). 

Integrating over t, letting p — oo, and recalling f = m^;, we have 

sup < C||9,Mo||oo, (2.10) 

0<t<To 

for sufficiently small Tq = To(||uo||oo, ll^lb, ||'9x%||oo)- 

Step 3: Control of ||(9xxii(^)||p for 2 < p < +oo. First take 2 <p < 
+00 and compute 

2 j \dxx'U^ dxxU ' ^ xxxxdx -\- 6 J" \OxxU\^ dxxUUxxxxdx 

- J \dxxu\^~'^dxxudxxx{{dxK * u)u)dx. 
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We further estimate by using integration by parts, 

J" \dxx'^\^ dxxU ■ (it ^ xxxxdx 

=3 J \dxxu\^~'^dxxU ■ u^dxxxxudx + 24: J \ dxxu\^~'^ dxxU ■ udxudxxxudx 

+ 18 y \dxxu\^udxxudx + 36 y Idxxul^idxU^dx 
= — 3{p — 1) y {dxxxuY\dxxu\^~'^u^dx + 18(1 ~ P) J \dxxu\^~'^dxxU ■ udxudxxxUdx 

+ 18 [ Idxxul^idxufdx. (2.11) 



Using Cauchy Schwartz inequahty: 



1 3 
ab < -a^ + -6^, 
- 6 2 ' 



we bound the second term in (12. lip as 
|18(1 — p) J \uxxf~'^dxxUudxudxxxudx\ 

<27(p-l) J \uxx\'^\dxu\'^dx + 3{p - 1) J \uxx\^~'^u^\dxxxu\'^dx. 

Hence, 

j \dxxuf~'^dxxU ■ {u^)xxxxdx <{27p - 9) J \dxxu\^idxufdx. 
Also it is obvious that 

e I IdxxUr'dxxUUxxxxdx 

= -ep j {dxxxuf\dxxu\^~'^dx < 0. 

Similarly we have 

\dxxu\^~'^dxxU ■ dxxx{{dxK * u)u)dx 

+ / \3.M'-'-a..rS,ua,K , a..ud. + / \a,,ur'a.,uua.j< , a„M. 
< j \a„u\'a.K , a^ud. + / \a„ur-a„ua.ua^K , a„ud. 

+ j \dxxu\^udx + j \dxxu\^~'^dxxUu{e"^^'^ * dxxu)dx 
<\\dxxu\\l-{\\dxu\\^ + \\u\U. 



< 
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Collecting all the estimates and we have 



\U\ 



Integrating over t we have 



sup ||9,xu(t)||p<e^P^"||9,,no||p. (2.12) 

0<t<To 

To estimate ||u^x||oo we need to estimate ||u^.^.^.||2. We have 



d_ 

dt 



r 

/,.4,H..„..:,.=K.-/.„„»,....K,* 



< 



Since 



< 



dxxxx{{dxK * u)u)uxxxdx 



^||^2;2;x||2ll^x^||oo + || "^rax || 2 || *9xa;^ || 2 || ^^xxx^ || 2 II ^ II oo- 



Finally we get 

d 



II2 ~ ||^x'3;3;||2(ll^x'||2 ||^x||oo + ||^||oo; 



-L II II llr) l|2 I II l|4 _L II l|2|| ||2 

H~ ll'^xxx II 2 II Cxx'^ll 2 ~r ll'^x'X'IU + || "^x-x || 2 II ^a; || 00 ' 

Gronwall then implies that 

sup ||9a;a;xM(i)||2 < C'||9:i.^^Mo||2- (2.13) 

o<t<r 

Since ||u^.x||oo < ||u||2 ||<9xxx^i||2% we conclude that 

sup ||Ua;^(t)|| 00 < C. 

0<t<To 



- / u" 
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Step 4: Control of ||9^'''"'"^m||2. We compute 

- j df+^{K * d^uu)df+\dx 
=: I + II + III 
For /, we integrate by parts and obtain 

J2 / diud^udydl^^\dx 

l+m+n=2k+2 
l<m<n<2k 

{df^'^ufdx + ^ j{ul + uu^^){df^^ufdx (2.14) 
+ Q,m,n I dlud^udx^'udf+^udx (2.15) 

l+m+n=2k+3 
l<m<n<2k+l 

Discarding the negative term and using Holder's inequality, we bound 
(12A411 as 

(EH < ilMlo + ||w||oo||Mx.x||oo)i|9f +^M||2- (2.16) 

To estimate the summand in (12.141) . we discuss several cases. 

Case 1. n = 2k + l. In view of the constraint, we have (/, m) = (0, 2) 
or (1, 1). Holder's inequality then gives that 



- / 



J diu&:udydf^'udx\ < (Ml + ||n|U||n,.,.|U)||9f +1 



Case 2. n < 2k, I = 0. Since l + m + n = 2k + 3 we must have m > 3. 
We choose 2 < n, g < oo such that - + - = ^ and use Holder to get 

p q z 

diud^udyaf+'udxl < ||n|U||9>||,||a>||J9f 
Using the interpolation inequality 



m - 2 - i 

p 



we have 

1 1 d[ud^ud:udf+^udx\ < \\u\U\df+M\lM^'i 

< Moo\\u\\HA\df^M\l 



14 



DONG LI AND XIAOYI ZHANG 



Case 3. n < 2k, 1 = 1. In this case m, n > 2, we can choose 

- + - = - 



2 < p,q < CO such that - + - = ^ and use the interpolation inequahty 



- 1 -^-^ 



||a>||, < ||n|| .3^^- (2.17) 

/3 2 



to get 



< Ik.llooll^f +^ 

< \K\U\u\\H4df^' 



u\\o\\u\\ ■ 3 
II2II Wfj^ 



M n. 



Case 4. n < 2k, l,m,n > 2. Choosing 2 < p,q,r < oo such that 
1+1+1=1 



^ + - + - = I and using the interpolation inequality fl2.17p . we get 



diud:^udydl'+\dx\ < ||ai«||,||a>||ja>|M|af +^ 



U\\2 



'^5 " ^'^^ 



Collecting all the estimates together, we conclude 

/ <(||9^m||^ + ||m||oo||mxx||oo + ||m||oo||m||h3 
+ ||mx||oo||m||h2 + lkll?/2)||9^''+^M||2- 

For term //, we simply have 



II = -e J {df+^ufdx < 0. 
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For III, we compute: 

III=- j df+\K * d^u)udf+\dx - j K* d^udf+\df+\dx 
+ Ci^^ j di{K * d^u)d^udf+^udx 



l+m=2k+2 
l,m<2k+l 



J d,,K * df+\udf+\dx + ^ j d,,K* d^Mdf^^ufdx 
J2 Ci,m I d^K * diud'^udf+^udx 

l+m=2k+2 
l,m<2k+l 

2 J u{df^\fdx - j (e"l"l * df^^u)df^\udx + i y {d^K * d^u){df^\fdx 
Ci^m J dj{ * diud'^udf+^udx 



l+m<2k+2 
l,m.<2k+l 



i2fe+l„.l|2 



<i\\u\\oo + \\dM2)\\di'^'u 



J2 Ci^mWdl'-^'uUdiuhWd: 



l+m=2k+2 
l,m<2k+l 



Using the interpolation inequality 
we finally get 

in<{\\u\\^ + \\u,h)\\dl'^'u\\l. 

Summarizing the estimates of /, //, /// we obtain 



+ \\u\\oo\\u\\h3 + kx- oo + n oo + Fx 2 



■i + \\U\ 



Using Gronwall and (EID, fl2A2ll . fl213l) we get 

sup \\df^'u{t)h<e^^^\\df^\42. 

0<t<To 

for sufficiently small Tq = To(||mo||2, ||<9xWo||2)- 

This concludes the estimate of the i7^^~'"^-norm of u"^. Finally if 
Mo > 0, then by the weak maximum principle we have u'^(t) > for 
any < t < Tq. The lemma is proved. □ 

We are now ready to complete the 
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Proof of Proposition \2.1[ This follows directly from Lemma 12.51 and 
12.61 In particular note that in Lemma 12.61 the time of existence of 
the local solution does not depend on e. □ 

Proof of CorollarylKM Assume Uq G CQ^f^H^{R)f]LP{R) for some 
1 < p < 2. Then using (12. ip and Duhamel's formula, we can write 

u^it) as 

- [ e'^'-'^^''^d^{u'd^K*u'){s)ds. 
Jo 

Since ||e^*^'"''/||p < for any e > 0, we can estimate the norm of 



h'{t)\\p<\\uo\\p- 

<||mo||p+ / \\u^{s)\\w^.3p + \\u'{s)\\'^i,2pds 
Jo 

fh'is)rm + \\u%s)\\i,ds, 

Jo 



ds 

p 



< ||mo|Ip + 



where the last inequality follows from the Sobolev embedding. This 
estimate shows that u'^{t) G L^. for any t. The continuity (including 
right continuity at t = 0) follows from similar estimates. We omit 
the details. Finally since for 1 < p < 2, ||uo||2 + ||'9a:Mo||oo ^ ll^ollp + 
||f^i:'Wo||oo; one can choose the time interval sufficiently small depending 
only on ||uo||p + ||9^Mo||oo- □ 

2.2. Proof of Proposition 12.31 In the case uq > 0, we shall show 
the regularized equation has a global solution. The key point here is 
that by using the positivity, we can obtain the apriori boundedness of 
the L^-norm {2 < p < oo) of the solution on any finite time interval. 
To control the L!^-norm of the gradient, we need the following lemma 
from [391. 



Lemma 2.7. Let Q = (0,T) x M and u G C2'i(fi) nC(fi) he a solution 
to the parabolic equation: 

dtu = b{t, X, n, n^O^xx + a{t, x, u, u^), {t, x) G (0, T) x M 

u\t=o = Uq 

with the following properties: 

(1) Uniform parabolicity : there exists constants /3o, (3i > such 
that 

l3oh(t,x,z,p)p^ > \a{t,x,z,p)\, V|p| > A. 
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(2) There is a constant M > such that 

\uit, x) - u{t, y)\<M, V (t, x), (t, y) G 

(3) There exists Lq > such that 

\uq{x) - UQ{y)\ < Lq\x - y\, Vx,?/G]R. 
Under all the above assumptions, we have the following hound: 
sup ||M,(t)||oo <2(Lo + /3i)e**'^ 

0<i<T 

Proof. See Lemma 11.16 of [M]- D 

Incorporating Lemma 12.71 with the idea of the proof of Lemma 12.61 
we get the following global analogue of Lemma 12.61 

Lemma 2.8. Let uq G H^''~^^{R), k > 2 and uq > 0. Then for any 
T > and for any e > 0, 02.11] has a unique solution u'^ G B^. Also 
u'^{t) > for any t > 0. 

Proof. The positivity of ti^ follows easily from the weak maximum prin- 
ciple. By Lemma 12.51 we only need to get an apriori control of the 
H^'^'^^-noim of u on an arbitrary time interval [0,T]. 

Step 1: L^-norm control for 2 < p < +oo. By the same estimates 
as in Step 1 of the proof of Lemma 12.61 we have the following a priori 
L^-norm estimate: 

\\u{t)\\2 < llnolhexp 

To obtain the L!^-norm estimate, we take 2 < p < +oo and compute 

-^\\u(tm<- '^^^'^^ I u\u?dx + ^^ I {e-\''\*u)\uV'dx 
pdt ^ p J p J 

Since m > we can drop the first term and continue to estimate 

~^\\umi<C\\u4,e^'\\um% 

where C is a generic constant independent of p. Gronwall then implies 
that 

\Ht)\\p < llwollpexp {Ct\\uo\\2e^') . (2.18) 

Letting p oo, we obtain 

IkWIloo < llMolUexp (CtllMoIhe^*) . (2.19) 

This is the L^-norm estimate we needed. 

Step 2: Control of ||92;M(t)||oo- We shall apply Lemma [2.71 with 

b{t,x,z,p) = 2;^ + e, 

a{t, X, z, p) = 2zp'^ + 2z^ — c{t, x)z — d{t, x)p, 
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where 

c{t,x) = {d^K * u^){t,x), 
d{t,x) = {dxK * u){t,x). 

By step 1 we have 

\c{t,x)\ + \d{t,x)\ < C\\u{t)\\oo 
< c\\ II oo 

exp{CT\\uo\\2e ). 

Also 

\u{t,x) - u{t,y)\ < 2||n(t)|U < 2||no||ooexp (CT||no||2e^^) . 
=: M 

and 

\uo{x) - uo{y)\ < ||5^-Mo||ook - y\ 
=: Lo\x - y\. 

Collecting all these estimates, we see that 

sup \\dxu(t)\\oo < C(||no||2, \\dxUo\\oo,e,T) 

0<t<T 

This concludes the gradient estimate. 

Step 3: Control of the higher derivatives. This part of the esti- 
mates is exactly the same as the corresponding estimates in the proof 
of Lemma 12.61 Note in particular that we do not need T to be small 
once we obtain the a priori control of ||m||p and ||5^m||oo norms. □ 

Now we are ready to complete 

Proof of Proposition \2.3[ This follows directly from the local existence 
Lemma 12.51 and the a priori estimate Lemma l2.8[ If uq G i^S(IR) for 
some 1 < p < 2, then by repeating the proof of Corollary 12. 2l we obtain 
u G C([0, T], LP) for any T > 0. The proof is finished. □ 



3. Proof of Theorem I1.1L 11.31 and 11.41 

In this section we prove our main theorems. Our solution is going 
to be the limit of the sequence of regularized solutions which we 
constructed in the previous section. To obtain uniform control of the 
Sobolev norms of the regularized solutions, we have the following 

Lemma 3.1 (A priori estimate). Assume uq G H^^^^iM), k > 2. Let 
To = To(||no||2 + ||<9x%||oo) o-nd u" G B!^^ be the corresponding unique 
solution to (12.11) (see Lemma W^) . Then the set of functions (M'')o<e<i 
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satisfies the following uniform estimates: 

sup \\u'\\H2k+i + Ve\\u'\\L2H2k+2nQTo]xm < C {\\Uo\\ H2k+i , k) , 
0<t<To 

||u''(9^''"^^U^||i2^2([o,To]xM) + \\dtU''\\L2H2k{[o,To]xR) < C {\\uo\\ H^k+i , k) , 
\\d^{d^K * u'u')\\L2H2k{[o,To]xM) < C{\\uo\\H2k+i,k), 

where C = C{\\uo\\jj2k+i,k) is a positive constant. 

Proof. The estimates of II II j;^ooj|^2fc+i, A/e||M''||^2j|^2fc+2 and ||m^<9|^"'"^m^||l2^2 
can be recovered from step 4 of the proof of Lemma 12. 6[ To bound 
||9tM''||2,2j:^2fe, we use fl2.1l) to obtain (here we again drop the superscript 
e for simphcity): 



\H2k < Wdtuh+Wdl'^dtuh 

< 11(^2 + e)u^^ + 2uul -d^{K* d^uu)\\2 + \\df{u\^,^)\\2 
+ e\\df+\h + ndlHuuDh + \\df^\K * d^uu)h 

< \\u\\%2 + ||m||^2 + ||M^9f +2^112 + ||M||H2fe+i 
+ e||9f +\||2 + ||M||H2fe+i + M\H2k+i- 

Integrating over [0,To], we have 

/•To 

/ ||c^i^i(^)|lH2fe'^^ ~ ^0 sup ||M||^2fe+i(l + ||n||//2fc+i) 
Jo 0<t<To 

+ r \\u^df+Mmldt + e r \\df^M\ldt 
Jo Jo 

< C{Tq, k, ||uo||i^2fc+i). 

Since Tq = To(||no||2 + llf^x^iolloo), we obtain the desired bounds. The 
last estimate is a simple application of Holder and Young's inequality 
using the given estimates. □ 

Now we define the set 

■. = [v e C([0,To];iJf ) nL°°(0,To;iJf+i)) : 

dtv e L\0,To;H^'')), v'd,,v e L'{0,To; H"')y 

We shall prove the existence and the uniqueness of the solution of (11. ip 
in this set. This is 

Theorem 3.2. Assume uo G H'^''^^{M.) with k > 2. Then there exists 
Tq = To(||no||2 + ||f?a;Mo||oo) such that (II. ip has a unique solution u G 

Proof. Step 1. Contraction in C([0, T], L^). For any two ei,e2 > 0, 
let u = u'^'^, V = m'^2 solve (12.10 with the same initial data Uq. Denote 
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w = u — V, then for w we have the equation 

dtw = ^{u^ - v^)^^ + eiu^^ - e^Vxx - dx{K * d^wu) 

Then 

— j dx{K * dxwu)wdx — j dx{K * dxvw)wdx 

=:I + II + III + IV. 

Denote F{u,v) = u"^ + v"^ + uv. Clearly we always have F{u,v) > 0. 
We compute 



/ = ^y dxx{wF{u,v))wdx 



dxxWwF^u, v)dx 



_ 1 
~ 3 

= y {dxwfF{u,v)dx+^ j w^dxxF{u,v)dx 

-\ j w'^dxxF{u,v)dx 

< \H\l{\\u\\l^ + \\v\\l^). 
For II we have the estimate 

II < {ei\\u\\H2 + e2||t'||//2)||u;||2. 
The term III can be bounded as: 

/// = [ K * dxWUWxdx 

dxxK * wuwdx — J K * dxWUxwdx 

< i\\u\\oc + ||5x'M||oo)||w||2. 

Similarly for IV we get 

IV<\\K*dxxV\\oo\H\l 
^ ll'y||Hl|kll2- 

Finally we have by Lemma 13.11 

+ (ei||u||^^2 +e2|K;||//2)||i/;(t)||2 
< C{\\uo\\H2,+i,k)\\w{t)\\l + C{\\uo\\H2k+i,k){ei + €2)\\w 
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Gronwall then gives 

\\u^^{t)-u^^{t)h = \\w{t)h 

< C{\\uo\\H2k+i,k){ei + 62) exp {C{\\uo\\H2k+i,k)t) . 

This imphes that (u^)o<e<i has a hmit as e ^ in C([0,To]; L^). De- 
note this hmit as u G C([0,To]; Ll). 

Step 2. We show that u is our desired solution in By Lemma 
13.11 and the interpolation inequality 

1 2k 



we obtain 



u'^u inC([0,To]; i/f ), as e 



It then follows easily that 

{u^fd^^u^ u^d^^u, in C(0,To; L^), 
u%d.Xf ^ u{d,u)\ in C(0, To; L^), 
d,,{K * d.Xu') ^d,{K* d^uu), in C(0, Tq; L^). 

By Lemma 13.11 the functions {u'^Yd^xU", u'^^d^u'')^, dx{K * d^u^u") 
are uniformly bounded in L^([0, Tq], ifj^). We conclude that u'^d^xU, 
u{dxu)'^, dx{K*dxUu) e T^([0, Tq],HI^). Now since u!^ is also uniformly 
bounded in L~([0, Tq], i/f +^), we obtain u e L°°([0, Tq], i/f Simi- 
larly by Lemma [3?T1 the set of functions (i9tu'^)o<e<i is uniformly bounded 
in T^([0, To], Hi''). By extracting a subsequence if necessary, we obtain 
for some e„ — > 0, 

d^u'^^^d^u, inL2([0,To],iff). 

It follows that dtu G T^([0, Tq], i^J'^). Finally by using Lemma [HT] again 
we have 



endxxu'-^Q, mL^{0,To;H 



2k\ 
X ) 



This shows that u is our desired solution in D^^. The theorem is 
proved. □ 

Corollary 3.3 (Continuation of solutions). Let uq G if^^"*"^ (M) with 
k > 2. Then there exists a unique solution u of (11. ip with maximal- 
lifespan T* such that only one of the following possibilities occur: 

(1) T* = 00 and u E D!^ for any T > 0. 

(2) T* < 00, ueD^ for any T < T* , and 

lim \\dxu(t)\\oo = 00. 
For any 2 < p < 00, there exists a generic constant C such that 
\Ht)\\p< hollpe^P*, Vt >0. 
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If in addition Uq > 0, then u{t) > and we also have the p-independent 
estimate for all 2 < p < +oo; 

\\uit)\\p<\\uo\\pexp{Ct\\uo\\2e^'), Vt > 0. (3.1) 

In particular if p = +oo, then 

\Ht)\\^ < llnolUexp {CtWuohe^') , Vt > 0. 

Proof. This follows directly from Theorem 13.21 and Lemma 12.61 (see 
(12. 5p for the growth estimate of L^-norm for any T > 0). Since u 
uniformly in C([0,T]; -ffj^), u{t) satisfies the same estimate (12. 5p . Sim- 
ilarly if Mq > 0, then the estimate (13.11) follows from the corresponding 
estimate for u"^ (see (I2.18P and (12.190 in the proof of Lemma [2.80 . Now 
that the norm of the constructed local solution remains finite for 
any T > 0, it can be continued as long as ||9a;M(t) ||oo does not blow 
up. □ 



We are now ready to complete the 



Proof of Theorem \l.l\ By Theorem 13. 2^ uq G f]^^Qi7™(M) implies 
that u E nr=2^To- Therefore u G C{[0,To],H^) for any m > 0. It 
also follows that dtu G nm=i -^^(0' ^o! H^). By using (11.11) . a standard 
bootstrap argument implies that u G C°°([0,To] x M). If uq > 0, then 
the positivity of u follows from Lemma and the uniform convergence 
of u*^ to u. The theorem is proved. □ 



Proof of Theorem M.'A This follows immediately from Theorem 1 1 . 1 1 and 
Corollary 13.31 □ 



Proof of Theorem L4_. This is only slightly different from the proof of 
Corollary [221 Assume uq G CC^f^H^(R) f]LP{R) for some 1 < p < 2. 
Then using (11.10 . we can bound 



Wit)\\p ^ \\uo\\p+ [ 
Jo 



ms 



d^{{d^K * u)u) 



+ \\u{s)\\wi,2pds 



ds 



^ \\M\p + 

^ll'«o||p+ / \\u{s)\\]jz + \\u{s)\\]j2ds 
Jo 



This shows that u{t) G for any t. The continuity (including right 
continuity at t = 0) follows from similar estimates. We omit the stan- 
dard details. Finally in the case Uq > 0, Mq G -Z^;J;(ffi), the L\ preserva- 
tion follows from direct integration. □ 
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We argue by contradiction. Let uq G C^(M), uq > and assume 
that the corresponding solution u is global. By Theorem 11.11 11.31 11.41 
u e C([0,T), if™) for any m > and the L^-norm of u is preserved. 
Our intuition of proving the blowup is based on the observation that as 
time goes on, the boundary of the solution (which in ID case consists 
of two points) will move face to face at a certain speed which has a 
lower bound independent of time. This clearly will lead to the collapse 
of the solution. To realize this intuition, we will carry out a detailed 
analysis on the characterstic lines of the solution which satisfy 



By standard ODE theory and our assumption that u is a smooth global 
solution, X{t) is well defined and smooth for all time. Moreover, we 
have the following lemma. 

Lemma 4.1 (Properties oi X{t,a)). The characterstic lines X{t,a) 
and the solution u{t,x) satisfy the following properties: 

(1) For any t > 0, X{t, ■) : M — * M a C°° diffeomorphism. 

(2) X(t, ■) is an order-preserving map in the sense that if ai < a^, 
then X(t, ai) < X{t, 02) for any t >0. 

(3) X(t,-) maps intervals to intervals. More precisely, for any ai < 
a2, denote 



X{t, ■){[ai, 02]) = {X{t, a) : ai < a < 02} • 

r/ien X(t, ■)([«!, 02]) = [^(^,tti],^(^,a2)]- 

(4) Ifuo{a) = 0, then u{X{t,a),t) = for any t>0. 

(5) If supp{uq{-)) C [—L,L] for some L> 0, then 

supp{u{t,-)) C [X{t,-L),X{t,L)] C [-L,L] 
for allt > 0. 

Proof. Property (1) is rather standard. Note in particular that 



is smooth and globally invertible due to the fact that u is smooth. 
Property (2) is also trivial. Assume not true, then by the intermediate 
value theorem, for some r > we must have X{t, ai) — X{t, 0^2) = 
which is a contradiction to property (1). Property (3) follows imme- 
diately from property (1) and (2). For property (4), one simply uses 
(11. ip and (14. ip . Writing u(t,a) = u{X{t,a),t), we have 



4. Proof of Theorem 11.51 




j^Xit.a) = {K*d,u){X{t,a),t), 
X{0,a) = a G M. 



(4.1) 




f^u{t, a) = {-dji * d^u + ud^^u + 2{d.^uY) u{t, a), 
m(0, a) = Uo{a). 
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On the RHS of the above equation, one can regard the coefficient of 
u{t,a) as given functions of (t,a). By direct integration, property (4) 
now follows immediately. For property (5) we take y < x(t, —L), then 
by property (1) and (2) we have y = X(t, a) for some a < —L. By 
property (4) since uo{a) = 0, we must have u{t,y) = u{t,X{t,a)) = 0. 
Similarly once can show that ii y > X{t,L) then u{t,y) = 0. This 
shows that supp{u{t, ■)) C [X{t, —L),X{t, L)]. Now using this fact and 
(14. ip . we compute 

^X{t, -L) = {d^K * u){t, X{t, -L)) 
at 

sign{X{t, -L) - y)e-^''^''-'^^-y^uit,y)dy 

sign{X{t, -L) - 2/)e-l^(*'-^)-J^ln(t, y)dy 

ly>X{t-L) 

> 0. (4.2) 

By a similar computation one can show that f^X{t, L) < 0. These facts 
together with property 2 and 3 above easily yield that [X{t, —L),X{t, L)] C 
[-L,L]. □ 

Using lemma mH we are able to finish the 



Proof of Theorem \1.5i We argue by contradiction. Let uq G 
uq > with supp{uo) C [—L, L] for some L > 0. Assume that the 
corresponding solution u is global. Define the characteristic X{t,a) 
according to (14. ip . By Lemma 14.11 we have 

supp{u{t,-)) c[X{t,-L),X{t,L)] C [-L,L], Vt > 0. (4.3) 

By this fact, we then compute using (14. 2p . 

^X{t, -L) = - [ sign{X{t, -L) - y)e"l^(*-^)-^ln(t, y)dy 

"I Jy>X(t~L) 

sign{X{t, -L) - y)e-^''^''-'^^-y\u{t,y)dy 

X{t-L)<y<X{t,L) 

^-\Xit,-L)-y\^^^^y^^y 

X{t-L)<y<X(t,L) 

> e"^^||uo||Li, 

where in the last step we used the fact that the norm of u is pre- 
served. This shows that X{t, —L) grows linearly with time which is 
contradiction to (14.31) . Thus we have shown that the solution u with 
Uq as initial data cannot exist globally. Finally (11.41) and (II. 5p are easy 
consequences of Theorem 11.11 11.31 and 11.41 The theorem is proved. □ 
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